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We investigate the effect of free electrons on the quality factor (Q) of a metallic nanomechanical 
resonator in the form of a thin elastic beam. The flexural and longitudinal modes of the beam are 
modeled using thin beam elasticity theory, and simple perturbation theory is used to calculate the 
rate at which an externally excited vibration mode decays due to its interaction with free electrons. 
We find that electron-phonon interaction significantly affects the Q of longitudinal modes, and may 
also be of significance to the damping of flexural modes in otherwise high-Q beams. The finite 
geometry of the beam is manifested in two important ways. Its finite length breaks translation 
invariance along the beam and introduces an imperfect momentum conservation law in place of 
the exact law. Its finite width imposes a quantization of the electronic states that introduces a 
temperature scale for which there exists a crossover from a high-temperature macroscopic regime, 
where electron-phonon damping behaves as if the electrons were in the bulk, to a low-temperature 
mesoscopic regime, where damping is dominated by just a few dissipation channels and exhibits 
sharp non-monotonic changes as parameters are varied. This suggests a novel scheme for probing 
the electronic spectrum of a nanoscale device by measuring the Q of its mechanical vibrations. 
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I. INTRODUCTION 

The design and fabrication of high-Q mechanical res- 
onators is an ongoing effort that has intensified with 
the advent of microelectromcchanical systems (MEMS) 
and even more with the recent progression toward na- 
noelectromechanical systems (NEMS). 1-3 One requires 
low-loss mechanical resonators for a host of nanotechno- 
logical applications, such as low phase- noise oscillators; 4 
highly sensitive mass, 5-8 spin, 9 and charge detectors; 10 
and ultra-sensitive thermometers 11 and displacement 
sensors; 12 ' 13 as well as for basic research in the meso- 
scopic physics of phonons, 1 1 and the general study of the 
behavior of mechanical degrees of freedom at the inter- 
face between the classical and the quantum worlds. 15-18 
It is therefore of great importance to understand the 
dominant damping mechanisms in small mechanical res- 
onators. 

A variety of different mechanisms — such as internal 
friction due to bulk or surface defects, 19-27 phonon- 
mediated damping, 2S ~ 3 1 and clamping losses 32-3 '' — may 
contribute to the dissipation of energy in mechanical res- 
onators, and thus impose limits on their quality fac- 
tors. The dissipated energy is transferred from a par- 
ticular mode of the resonator, which is driven exter- 
nally, to energy reservoirs formed by all other degrees 
of freedom of the system. Here, we focus our attention 
on electron-phonon damping, arising from energy trans- 
fer between the driven mode of the resonator and free 
electrons. This dissipation mechanism is avoided alto- 
gether by fabricating resonators from dielectric materi- 
als, but for different practical reasons one often prefers 
to fabricate MEMS and NEMS resonators from metals, 
such as platinum,' 5 ' gold, 38 ' 39 and aluminum. 40 Free 
electrons are also present in metallic carbon-nanotubc 
resonators 43 ' 44 and in resonating nanoparticles. 45-4 ' All 



these different resonators exhibit a wide range of quality 
factors, from as low as about 10 and up to around 10 5 , 
yet one still lacks a full understanding of their damping 
mechanisms. 

It is well-known from at least as early as the 1950s 
that electron-phonon scattering is a dominant source 
of attenuation of longitudinal sound waves in bulk 
metals at low temperatures, 48-53 and indications ex- 
ist that it may play a significant role in the damp- 
ing of longitudinal vibrations in freely suspended bi- 
pyramid gold nanoparticles. 4 '' We note that the ef- 
fect of electron-phonon scattering on electronic trans- 
port through suspended nanomechanical beams, 51 car- 
bon nanotubes, fullerenes, 58 atomic wires, and 
molecular junctions '" - '''' is well documented and in- 
tensively studied. There is also evidence for the ef- 
fect of electron-phonon scattering on heat transport in 
nanostructures.'"' 08 Motivated by all of these consider- 
ations, it is our aim here to estimate the contribution 
of electron-phonon interaction to the damping of vibra- 
tional modes in small metallic resonators, while focusing 
on the effects of their finite dimensions. 

We describe the interaction between electrons and 
phonons by means of a simple screened electrostatic po- 
tential. We assume that initially both electrons and 
phonons are at thermal equilibrium at the same tem- 
perature, except for a single mode, which is externally 
excited by the addition of just a single phonon to its ther- 
mal population. This allows us to assume that the elec- 
trons remain almost thermally distributed at all times, 
even though they do not actually relax back to equilib- 
rium. The decay rate of the excited mode is calculated 
pertubatively, using Fermi's Golden Rule, as the differ- 
ence between the rates at which phonons enter and leave 
the excited mode through their interaction with free elec- 
trons. This requires us to assume that the electron and 
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phonon energies are precisely known, or in other words 
that the a priori lifetimes of both the electrons and the 
phonons are much longer than all other relevant time 
scales. 

For the phonons this means that all other damping 
mechanisms must be much weaker than electron-phonon 
damping — although we show later that additional damp- 
ing mechanisms do not significantly alter the results of 
our calculations. For the electrons, on the other hand, 
this implies that we are working in the high-frequency, 
or unrelaxed adiabatic limit, with w q T e > 1, where uj q is 
the vibration frequency and r e is the mean lifetime of the 
electron due to its scattering with other electrons, ther- 
mal phonons, defects, etc. In certain situations, as dis- 
cussed in detail in section 5.12 of the book by Ziman, J ~ it 
is sufficient to satisfy the spatial version of this require- 
ment, namely that A e q > 1 — where A e is the electron 
mean free path and q is the wavenumber of the excited 
mode — which is easier to satisfy because the ratio of the 
phonon group velocity to the electron Fermi velocity is 
typically very small. Intuitively speaking, it is as if the 
moving electron explores the elastic wave much faster 
than it would have, if it were standing in place and wait- 
ing for the wave to go by. In bulk metals it is difficult 
to satisfy the adiabatic condition, and one is inevitably 
required to address the relaxation of the electronic dis- 
tribution by employing the Boltzmann equation or other 
approaches. 50 ' 51 ' 69 However, in clean nanometer scale de- 
vices oscillating at very-high frequencies, and operating 
at sufficiently low temperatures, there is a greater chance 
of reaching the adiabatic limit. We therefore assume that 
this is the case, and alert the reader to the fact that our 
results may be less applicable at high temperatures. 

We describe the vibrational modes of the nanomechan- 
ical resonator using continuum elasticity theory, which 
is often employed for treating nanomcchanical systems, 
even in the case of carbon nanotubes, l0 ~ 72 , and also in 
the quantum regime. 73-1 5 The small size of such nanos- 
tructures may raise the question of the validity of a con- 
tinuum elastic approach. However, explicit comparisons 
with atomistic calculations and experimental results have 
shown that continuum elasticity is valid down to dimen- 
sions of a few nanometers/ and may indeed be used 
even for carbon nanotubes, as long as one uses appropri- 
ate effective parameters. 80 ' 81 

Finally, we investigate beams with typical dimensions 
that are much larger than the bulk Fermi wavelength. In 
such systems it is usually assumed that the effect of the 
boundaries on the free electrons can be ignored and that 
the electrons are essentially unconfined. However, as we 
demonstrate here, the transverse dimensions of the beam 
set a temperature scale — between a few kelvins to more 
than a hundred kelvins for the beam geometries consid- 
ered here — below which the confinement of the electrons 
can no longer be ignored, and electron-phonon damping 
is expected to behave qualitatively different. 

The paper is organized as follows. In section II we 
describe the quantized vibrational modes of a nanomc- 



chanical beam. The interaction between these modes and 
free electrons is addressed in section III, and the result- 
ing expressions for the dimcnsionlcss damping Q _1 arc 
derived in section IV. The reader who is not interested in 
the technical details of the calculations is invited to skip 
directly to section V, where we discuss and explain the 
results of the calculation and their physical implications. 
We conclude with a summary in section VI. 



II. QUANTIZED VIBRATIONAL MODES OF A 
NANOMECHANICAL BEAM 

We consider a mechanical resonator in the form of a 
thin elastic beam with a rectangular cross section, al- 
though we expect the essence of our calculation to be 
independent of the specific geometry of the resonator. 
In what follows we describe the classical longitudinal 
and flcxural normal modes of vibration, obtained within 
standard theories of thin elastic beams, ' and quantize 
these modes as a collection of simple harmonic oscillators. 
We do not consider the twist modes of the thin beam be- 
cause the displacement fields of these modes have zero 
divergence and therefore do not couple significantly to 
free electrons, as explained in the next section. 

We take the length I of the beam to lie along the x axis 
from x = — 1/2 to x — 1/2, and its transverse dimensions 
along the y axis from y = to y = a, and along the z 
axis from z = to z = b. In equilibrium, the beam is un- 
strained, unstressed, and at temperature T everywhere. 
Departure of the beam from equilibrium is described by a 
displacement field lit (i = x, y, z). The displacement field 
Ui , as well as the strain and stress tensors Uij and are 
all functions of position and time. Yet, we assume that 
the temperature of the beam remains uniform and con- 
stant during the vibration, thus ignoring thermoelastic 
effects. 2<s We take the surfaces of the beam to be stress 
free, which implies that all but the a xx component of the 
stress tensor vanish on the surface. Because the beam 
is thin, this approximately holds in its interior as well. 
Hookc's law for the thin beam then takes a rather simple 
form: 

1 

^xx — ~f^O~xx: (1^7 

hi 

a 

Uyy = U zz = - —<Jxx, (lb) 

hi 

u xy = u yz = u zx = 0, (lc) 
where E is Young's modulus, and a is Poisson's ratio. 



A. Quantized longitudinal modes 

To describe longitudinal vibrations in the thin beam we 
make an additional simplifying assumption of neglecting 
any contribution to the dynamics that arises from having 
a nonzero Poisson ratio. We therefore take u y = u z = 0, 



3 



take u x to be independent of y and z, and ignore any devi- 
ation of the cross section of the beam from its rectangular 
shape. Under these assumptions, longitudinal vibrations 
are governed by a standard dispersionless wave equation. 
Thus, a mode of wavenumber q vibrates at a frequency 
uj q = ciq, where ci = y/E/p is the speed of sound of lon- 
gitudinal waves in the bulk, and p is the mass density of 
the beam. Taking a 7^ in Eq. (lb) would lead to the 
so-called Love equation 82 and to dispersive longitudinal 
modes, which we do not consider here. 

In the limit I — > 00 of an infinitely long thin beam 
it is convenient to describe the longitudinal modes as 
travelling waves leading to a phonon field operator of the 
form 



E 



2Mw„ 



(b g + bl q ), 



(2) 



where b q and b q are bosonic creation and annihilation 
operators, M = pV is the mass of the beam, V = abl 
is its volume, and we use a large volume normalization 
of the phononic wave functions as M and V tend to in- 
finity along with I. In this limit, the sum ^2 q is to be 
interpreted as an integral / j dq/2w. 

For beams of finite length we consider either doubly- 
clamped boundary conditions, with u x = at both ends 
of the beam, which is a common experimental geometry, 
or stress- free boundary conditions, with du x /dx = at 
both ends, which may be suitable for describing freely 
suspended resonators in solution. 45-47 The phonon field 
operator is then given by 



E 



V M ^q, 
n— 1 V H7 



(A n shvq n x + B n cosq n x) (b qn + b q J 



where q n = mr/l. and for doubly-clamped modes 



[A n ,B n ] = <J 

while for stress-free modes 
[A n ,B n ] = J 



[0,1], odd n, 
[1,0], even n, 



[1,0], odd n, 
[0,1], even n 



(3) 



(4) 



(5) 



rectangular shape of the cross section remains unaltered 
during the bending motion. Such an approximation is 
justified for small deflections since the error it introduces 
is only on the order of the transverse beam dimension di- 
vided by the radius of curvature of the bending. Since we 
assume a non-dcformable cross-section, there is in fact a 
neutral surface running through the length of the beam, 
at z — b/2, which suffers no extension or contraction 
during its bending. One can then show 8 '' that the lon- 
gitudinal strain component u xx , a distance Sz = z — b/2 
away from the neutral surface, is equal to Sz/R. By re- 
placing the curvature of the beam l/R with —d 2 w/dx 2 
we express the non-zero components of the strain field in 
the beam as 



d 2 w 
dx 2 ' 
b 

~ 2 



d 2 w 
dx 2 



(6a) 



(6b) 



Using this strain field it is possible to write down the 
Lagrangian density of the beam and derive its equation 
of motion. This equation contains two kinetic terms, 
one that is associated with the transverse motion of the 
beam and one that is associated with the rotation of the 
cross-section. The latter is smaller by a factor of b/R, 
and is therefore usually neglected, leading to the Eulcr- 
Bernoulli equation of motion 



pA 



d 2 w 



dx 2 



EI 



r d 2 w 
dx 2 



= 0, 



(7) 



where A — ab is the area of the cross-section, and / = 
a6 3 /12 is the moment of inertia of the cross section. The 
resulting flexural modes possess a quadratic dispersion, 
w q = Gq 2 , where G = yjEI/pA = b^E/Up. 

In the limit of an infinitely long beam, as in Eq. (2), it 
is again convenient to describe the transverse modes as 
travelling waves with a quantized phonon field operator 
of the form 



w 



E 



2Mw a 



f qx (b q 



6 f 



(8) 



For beams of finite length we consider doubly clamped 
boundary conditions, taking w = dw/dx = at both 
ends. The phonon field operator is then given by 



B. Quantized flexural modes 

We describe flexural vibrations in the thin beam by the 
transverse motion w(x, t) of the neutral axis of the beam 
in the z direction. We make the usual Euler-Bernoulli 
assumption that the transverse dimensions of the beam, 
a and 6, are sufficiently small compared with the length 
I of the beam and the radius of curvature R of the bend- 
ing that any plane cross section, initially perpendicular 
to the axis of the beam, remains perpendicular to the 
neutral axis during bending. We further assume that the 



w 



E 



„ 2Muj a 

n—l V H1 



(A n sinq n x + B n cosq n x 



+ C„ sinh q n x + D„ cosh q n x) (b qn + b qn ) , (9) 

where the wavenumbers q n are solutions of the transcen- 
dental equation cos q n l cosh q n I = 1, with q n l tending to 
odd- multiples of it/2 as n increases. The numerical coeffi- 
cients A n ,B n , C n , D n are determined through the bound- 
ary conditions and the normalization of the phonon wave 
functions, where by symmetry A n = C n — for odd n, 
and B n = D n = for even n. 
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III. ELECTRON-PHONON INTERACTION 

We assume a standard screened static interaction 
potential 84 



v (r — r') = 



-inZe' 

2 

<7tf 



S(r-r') 



(10) 



between the negative charge — e of the electron density 
So- V4-( r )Vv( r ) and the positive charge Ze of the distur- 
bance in the density no5v(r') of the ionic background, 
induced by the vibration. Here q TF is the Thomas-Fermi 
wavenumber, e is the magnitude of the electron charge, 
a is the electron spin (not to be confused with Poisson's 
ratio), Z is the number of valence electrons per atom 
in the material, no is the atomic density, and 8v(r) is 
the local volume change induced by the vibration. The 
electron-phonon interaction Hamiltonian, derived from 
this potential, is 



- P h 



r f^^(r^(r)fe(r),f 3 r, (11) 



where 



AnZn e 2 



(12) 



For small deformations, to first order in the displace- 
ment field, the local change in volume Sv(r) = V ■ u. 
For twist modes V ■ u = giving 5v(r) = 0, which is 
why we have ignored them here; for longitudinal modes 
V • u = du x /dx, where u x is given by either Eq. (2) or 
Eq. (3); and for flexural modes we find from Eqs. (6) that 



Vu = (2(7- 1) 



d 2 w 
dx 2 ' 



(13) 



where if is given by either Eq. (8) or Eq. (9). 

We take two approaches for describing the electron 
field ip(j), where from here onwards we suppress the spin 
index. As discussed in section I, treating the electrons 
as bulk-like is a common approximation for structures 



of the size considered in this work. Thus we can use a 
simple free electron field of the form 



m = J2 



i 



ik-r 

e c k , 



(14) 



where Ck is a fermionic annihilation operator. Within this 
approximation it is permissablc to treat the transverse di- 
mensions of the beam as being infinite from the point of 
view of the electrons. Accordingly, the sum is re- 
placed by a 3-dimcnsional integral V f dk/ (2ir) 3 . This 
leads to a simplification of the calculations and the re- 
sulting expressions below. 

Alternatively, we do not neglect the finite transverse 
dimensions (but the length of the beam is still considered 
to be infinite as far as the electrons are considered), and 
view the electrons as being geometrically confined. We 
then obtain the "particle in a box" field operators 



tp(x, y, z) = ^_ e «fc*z gm gm ^ z Ck ^ 



(15) 



with k y = riyTr/a, k z = n z Tt/b, and n y ,r 
both cases the energy of the free electrons is given by 



1,2,3, 



In 



2J.2 



2m 



K 2 k 
2m 



(16) 



but when the lateral confinement of the electrons is not 
neglected it takes the explicit form of parabolic bands 



tfkj 
2m 



■, with 




(17) 



A. Interaction of free electrons with longitudinal 
phonons 

To obtain the interaction Hamiltonian (11) for lon- 
gitudinal phonons with unconfined free electrons in an 
infinite beam we use the derivative of Eq. (2) for the 
quantized elastic displacement field, and the field opera- 
tor (14) for unconfined electrons. This yields 



oo 



tfe-ph = -*f £ J^q J dye< k »- k »> j dz^~ k > J d X ei k ^ + ")%c k (b q + bl q 

k,k',g " 



= -^EJ^7^ +q Ckk + ^), (18) 



k,g V 1 



where q = qx, and where we have performed the sum over exponential integrals in the first line of Eq. (18) . We note 
k' to eliminate the three delta functions that appear as that the integrations in the y and z directions involve only 
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the unconfincd electron wavenumbers. Thus, the width 
and thickness of the beam can effectively be taken to be 
infinite, which leads to momentum conservation in the y 
and z directions. 

When the confinement of the electrons is not neglected, 
the electronic field (14) is replaced with the one given 
by Eq. (15). When longitudinal phonons in an infinite 
beam are considered, the replacement of the continuous 
electron spectrum (16) with the set of parabolic bands 

I 



given by Eq. (17) yields the same formal expression for 
the interaction Hamiltonian as above in Eq. (18). The 
only difference is that the summations over k y and k z 
become discrete. 



To obtain the interaction Hamiltonian of longitudi- 
nal modes with unconfined electrons in a beam of finite 
length we use the derivative of the clastic displacement 
field given by Eq. (3). This then yields 



H, 



e— ph 



- y 

2.f^ V Mci 



dye 



i(ky-k' y )i 



dze 



2 1 k z —k„ 



dxe 



k x -k„ 



fA . D , . (q n + Ak x )l . (q n -Ak x )l 
(A n — iB n ) sine h (A n + iB n ) sine 



(-4 n cosq n x B n siuq n x) cj^Ck (pq n ~\~ ^J n ) 
c k+A k Ck(& 9 „+6tJ, (19) 



where Ak = Ak x x, with Ak x = k' x — k x , and where the 
sinc(a) = sin(a)/a functions replace the momentum con- 
servation delta function that exists for the infinite beam 
in Eq. (18). 

We note that by taking the limit of I — > oo in the 
first line of Eq. (19) we obtain the interaction Hamilto- 
nian for an infinitely long beam as written in the basis 
of standing waves instead of the basis of traveling waves, 
used in Eq. (18). Using this Hamiltonian it is possible 
to calculate the decay rate of a standing wave mode in 
an infinitely long beam, which is identical to the one ob- 
tained from the Hamiltonian in Eq. (18) for a traveling 
wave. Alternatively, by converting the creation and anni- 
hilation operators of the standing waves in Eq. (19) into 
creation and annihilation operators of traveling waves, it 



is possible to show that the Hamiltonian in Eq. (19) is 
identical to that of Eq. (18) in the limit of an infinite 
beam and continuous phonon wavenumbers. 

B. Interaction of free electrons with flexural 
phonons 

For flexural vibrations interacting with laterally con- 
fined electrons in an infinite beam, we substitute the di- 
vergence (13) of the quantized elastic displacement field 
(8) together with the field operator (15) for laterally con- 
fined electrons into the general expression (11) of the 
electron-phonon interaction Hamiltonian. We then ob- 
tain 



, a b oo 

ffe-ph = 4g(2 ^~ ^ J2 \/^![g J d V sink yy s ' mk 'yy J dz (z-^) sin k z z sin k' z z J dxe i ( k '-< + *)*<$ tf c k (b q + bl. 

k - k '.9 -oo 

_ 8g(l-2<r)b v / V n z n' z t / t \ 

where now Ak = qx + Ak z z, with Ak z = k' z — k z , and where the last sum is restricted to pairs of integers n z and n' z 
that are of different parity. Note that owing to the transverse nature of the vibration the electron momentum changes 
both in the x and in the z directions. 

The interaction Hamiltonian between the flexural modes and unconfined free electrons in an infinitely long beam is 

, oo oo b 

H c . ph = ^1 £ /Jg: / dxe^-^ J dye^~K)y J dz L - *) e^^U {h + ^) 



— oo 

g(l-2a) ^ / V 2sin^- Ak z bcos^ f / x \ 



, 2MG Ak 



where as above Ak = qx + Ak z z, with Ak z — k' z — the neutral axis of the beam, which is a result of the 
k z . The linear dependence of Sv on the distance from 
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Eulcr-Bcrnoulli thin-beam approximation, precludes the 
replacement of the finite limits of integration in the z 
direction with integration over an infinite range, even 
though we consider unconfined electronic states, since it 
leads to an unphysical divergence of the decay rate. A 
more realistic model for the flexural modes may saturate 



where here Ak 



replaces the momentum conservation delta function that 
exists in the infinite beam (20). 

IV. INVERSE QUALITY FACTOR 

The damping of mechanical vibrations is estimated by 
assuming that both electrons and phonons are in thermal 
equilibrium except for a single vibration mode, which 
is externally excited by adding just a single phonon to 
its thermal population. The rate at which the excited 
mode decays is evaluated as the difference between the 
rates at which phonons leave it and enter it due to their 
scattering with electrons. These rates are calculated with 
Fermi's golden-rule using the interaction Hamiltonians 
H c -ph derived above. The decay rate of an externally 
excited vibration mode of wavenumber q is thus given by 

r ? = y E \(.f-\H c _ ph \t)\ 2 S(e k -e k , + hio q ) 

11 k,k',<T 

~Y £ l(/+|tfc-ph||*)| 2 <5(£ k -£k<-^ 9 ), (24) 

k,k',(T 

where |/_) is the state with one less phonon and |/ + ) is 
the state with one more phonon, namely 

\i) = \n k ,n v ,N 9 + l), (25a) 
|/_> = |nk-l,n k / + l,JV,), (25b) 
|/ + ) = |n k -l,n k / + l,JV, + 2) > (25c) 

where, without any loss of generality, for infinitely long 
beams we consider modes with a positive wavenumber 
q > 0. The dimensionless damping, or inverse quality 
factor, of a mode with wavenumber q is then defined as 
Qq 1 = T q /uj q . In Eqs. (25) n k stands for the Fermi-Dirac 
distribution and N q for the Bose-Einstein distribution. 
We note that we do not take into account the change in 



this effect and lead to finite results even in the limit of 
large b. 

The interaction Hamiltonian for flexural vibrations in 
a finite beam is obtained by replacing the quantized dis- 
placement field of the infinite beam (8) with the one given 
by Eq. (9), yielding 



(22) 



(23) 



I 

the electronic chemical potential as the temperature is 
increased above zero, but rather assume that it remains 
equal to the Fermi energy of the beam. The latter is 
calculated for each specific thickness and width of the 
beam when the confinement of the electrons is taken into 
account. 



This procedure presupposes that the electron and 
phonon energies that appear in the delta functions in 
Eq. (24) are exact, or equivalcntly that the a priori life- 
times of electrons and phonons are infinite. For this to 
be valid we require (1) that all other vibration damp- 
ing mechanisms be much weaker than electron-phonon 
damping; and (2) that scattering rates of electrons — 
with other electrons, thermal phonons, defects, or surface 
imperfections — be much slower than the frequency of the 
mechanical vibration. As we demonstrate in section VD, 
removal of the first requirement does not significantly af- 
fect the contribution of electron-phonon scattering to the 
overall damping, as long as the different damping mecha- 
nisms are assumed to be independent of each other. The 
second requirement implies that our calculation is valid 
in the adiabatic limit, ui q T e > 1, where r e is the mean life- 
time of the electron. This requirement will be better sat- 
isfied for cleaner and smaller — hence, higher-frequency — 
devices at sufficiently low temperatures. 



In what follows we give the derivation of the exact 
expressions for the inverse quality factors, using Eq. (24) 
along with the interaction Hamiltonians of section III. 
The reader who is not interested in these rather technical 
derivations is welcome to skip to the next section where 
we discuss the physical consequences of these expressions. 



= Ak x x + Ak z z, and 
j_ 

2 

Ai ( K ) — J J elKX (An sm Qn x + B n cos q n x — C n sinh q n x — D n cosh q n x) dx 
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A. Damping of longitudinal vibrations by 
unconfined electrons in an infinite beam 

Using the Hamiltonian of Eq. (18) and the general ex- 
pression (24) for the decay rate we obtain 



Eq. (30) is independent of the mode wavenumber, the 
temperature, and the geometry of the beam. In fact, 
Q depends only on the material parameters of the beam, 
and typically varies between 10 and 10 3 for different met- 
als. 



Qa 1 



Kg 

Mcf 



k.cr 



[n-k (1 - ^k+ q ) {N q + 1)5 (e k - £k+ q + t^ q ) 



- n k (1 - rtk-q) (N q + 2)8 (e k - e k _ q - frw q )} , 

(26) 

where q = gx. We convert the sum over wave vectors in 
Eq. (26) into three integrals and change the energy con- 
servation delta functions into momentum delta functions 

<S( £k - £k±q ±^) = ^(^-(^T|)). (27) 

After performing the integration over k x , k v and k z and 
summing over the spin index we obtain 



2 2 



2tt h 4 pcfl3q 



1 



(N q + iy l^-ln- 



+ eT+ 



(N q + 2)N q I ptvjjq- In- 



where 



7± 



2m 



' cjjn q 

, h 2 



+ e'<- 



(28) 



(29) 



Both logarithmic terms in Eq. (28) are negligible as 
long as cim/H and q/2 are small compared to hp (which 
is the case for all reasonable values of q). For such vi- 
brational wavenumbcrs 7+ ~ 7_ and the logarithms van- 
ish. Combining the two terms in the square brackets in 
Eq. (28) then yields 



2 2 
g^m z 

2-nW'pci ' 



(30) 



which is the same as the result for the damping of 
bulk longitudinal acoustic waves in the adiabatic limit 
(Ziman, ,J Eq. 8.10.9, and Kokkcdce " '). The damping in 



B. Damping of longitudinal vibrations by 
unconfined electrons in a finite beam 



We substitute the Hamiltonian of Eq. (19) into Eq. (24) 
and again convert the sums over the electronic wave vec- 
tors into integrals. After summing over the spin index 
we obtain 



on f 

2 2 / dkdfc> k (1 - n k+Ak ) W„ (k x - k' x ) 
TT z pcfq n J 

(31) 



167T 2 /9Cj 

x [(N qn + 1) S (e k - e k +Ak + fiw g J 
- (N 9n +2)6 (e k - e k+A k - ^Jg n )] , 



where Ak = (k' x — k x )ic, and 

1A> ( \ A ( ■ ( K ~<ln)l . . (n + q n )l 
W n («) = A n sine h sine 



, . (K-q n )l . (n + q n )l\ 
B n I smc sine ) . (32) 



We perform the integration over k y and k z and change 
the energy delta functions into a delta function of the 
variable k' 



S (e k - Ek+Ak + fruq n ) 



H 2 k' x ° 

and a delta function of the variable k 
S (ek - e k +Ak - fwq n ) 



[6 (k' x -k' x °)+ S (k' x + k' x )}, (33e 



h 2 k° 



[S (k x -k° x )+S (k x + k° x )] , (33b) 



where (k' x ) = k x + 2mciq n /h, and (k x ) = k' x + 
2mciq n I h. We integrate over k' x for the first delta func- 
tion and over k x for the second delta function and obtain 



Qql 



2 2 
ng m 

8irh 4 pcfl3q„ 



(N q „ + 1) 

N qn (Nq n + 2) 



U/0 

h 



(3hu) qn + In 



J./0 



Pfmjq n + In 



M±( k 2_ k 2\ 
1 Jf. g 2r» { K ® K F) 

i + e 4£(*l+^-*S) 



1 



l + e 



: (*l~kl) 



[Wn (k x - k'?) + Wn (k x + k' x °)] 

[W n (k' x °~k x )+Wn(~k'°-k x )] 



(34) 
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where we have dropped the prime from the integration 
variable. 

The integrands in Eq. (34) nearly vanish for most val- 
ues of k x unless it is relatively close to ± (mci/h — q n /2). 
For these values of k x (and as long as q n is small com- 
pared to kp) both k 2 — kp and k 2 + 2mciq n /h — k\ are 
nearly equal to —k 2 , and the logarithmic functions in 
Eq. (34) can be neglected. Finally, by using the fact 
that W„(k) = W„(-k) we find that 

oo 
— oo 

+ W n (k x + k' x )] . (35) 

We note that the damping for the infinite beam in 
Eq. (30) and for the finite beam in Eq. (35), due to inter- 
action with unconfined electrons, are both independent 
of temperature and of the length of the beam for a given 
mode. 



confined electrons together with the general expression 
for the decay rate (24) yield 

q G2b2M k fe. - k ?f 

x [n k (1 - n k +q+Ak) (N q + 1) 

X 5 (e k - £ k +q+Ak + hug) 

- n k (1 - n k _ q+A k) (N q + 2) 

x 5 (e k - £ k - q +Ak - Twq)] , (38) 

where here Ak = (k' z — k z ) z. 

We sum over the spin index and change the delta func- 
tions in Eq. (38) into delta functions of the variable k x 

Tfl 

5 (e k - e k +Ak± q ± luJq) = T^-S (k x - K±) , (39) 



Damping of longitudinal vibrations by laterally 
confined electrons in an infinite beam 



As noted in section III, the Hamiltonian for the inter- 
action between longitudinal phonons and laterally con- 
fined electrons is formally identical to the one describing 
the interaction between longitudinal phonons and uncon- 
fined free electrons, with the integrals over k y and k z in 
Eq. (18) replaced by discrete sums over the allowed values 
These discrete sums yield 



of k y and k 



q; 1 9 m 



h 2 EAq 



n {ki,k y ^)^-^,k y ,k;))(N q + 2)] 1 (36) 



where 



h 2 



(37) 



D. 



Damping of flexural vibrations by laterally 
confined electrons in an infinite beam 



The Hamiltonian of Eq. (20) which describes the inter- 
action between flexural modes in an infinite beam with 



where 



K± 



mGq k\ — k' z 2 q 



± 



2q 



(40) 



We express the sum over k x as an integral, and use 
Eqs. (39) and (40) to obtain 



! _ 768. 9 2 to(2ct - l) 2 
' q ~ n 4 h 2 EAq 



n z n z 
^ [n 2 z - n' 2 ) 4 



x [n( K+ , kv ,K) (1 - n { K + + q ,k v ,k> z )) (N q + 1) 

- n {K _ ikyM {l-n {K __ q ^ kyM) ) (N q + 2)] .(41) 



E. Damping of flexural vibrations by laterally 
confined electrons in a finite beam 



The inverse quality factors of flexural modes in a finite 
beam are obtained in a similar way but without exact 
momentum conservation. Thus, the momentum delta 
functions are replaced with the C n functions, given by 
Eq. (23), and the integral over k x remains in the final 
expression 
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2m 
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2mu 
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2m 
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(42) 



£ 85 


70 GPa 


a 85 


0.35 


P 85 


2.7 gr/cm 3 




3 


\ 86 

Atf 


0.049 nm 



TABLE I. Material properties of bulk aluminum at room tem- 
perature, used for calculating the results in section V. 



V. RESULTS AND DISCUSSION 



We use the material properties of aluminum in our 
numerical calculations below, because, according to our 
model, it is expected to have the highest dissipation 
among the typical materials used in the fabrication of 
metallic nanomechanical beams. We use the room tem- 
perature values of these properties, listed in Table I, 
as they do not change much between room-temperature 
and cryogenic temperatures. ,V] We consider beams with 
cross-sectional dimensions that vary between lOnm and 
50nm, and finite lengths that vary between 500nm and 
2.0/im. Because our model involves many approxima- 
tions, the calculated values of Q^ 1 given below should 
be considered only as a rough order of magnitude esti- 
mate of electron-phonon damping. Nevertheless, we be- 
lieve that we describe correctly the qualitative behavior 
of the damping as a function of the different physical and 
geometrical parameters. 
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FIG. 1. (color online) Inverse quality factors, on a logarith- 
mic scale, of the first 50 longitudinal modes interacting with 
unconfined electrons, in a finite beam with stress-free and 
clamped boundary conditions, compared with Q~ l of an in- 
finite beam. 



A. Damping of longitudinal vibrations as a 
function of temperature and wavenumber 

The inverse quality factors of the first 50 longitudinal 
modes of a 1/im-long beam — calculated using Eq. (35), 
which assumes that the electrons are unconfined by the 
lateral dimensions of the beam — are shown in Fig. 1, for 
both clamped and stress-free boundary conditions. As 
one should expect, for sufficiently short wavelengths the 
damping is not affected by the finite length of the beam. 
Accordingly, both sets of results converge to the value of 
Q , calculated for an infinite beam using Eq. (30), which 
gives an estimated Q of about 40 for short-wavelength 
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FIG. 2. (color online) Inverse quality factor of the longitu- 
dinal mode q — 7r/10 3 nm _1 as a function of temperature in 
infinitely long beams with square cross sections, where the lat- 
eral confinement of the electrons is taken into account. The 
value of Q _1 for unconfined electrons is shown for comparison. 



(large q) modes in aluminum beams. For long-wavelength 
(small q) modes the finite length of the beam has a sig- 
nificant effect on the damping. The inverse quality factor 
of the lowest clamped mode is about 25% smaller than 
Q^ 1 of an infinite beam, and that of the lowest stress- 
free mode is almost 4 times larger than Q _1 of an infinite 
beam, with a difference larger than 10% persisting up to 
about the tenth stress-free mode. 

The lack of an intrinsic energy scale in the spectrum 
of the unconfined electrons — other than the Fermi en- 
ergy which is too high to be relevant — results in electron- 
phonon damping that is temperature independent for 
both infinite and finite beams, as is apparent from 
Eqs. (30) and (35). This changes when the lateral con- 
finement of the electrons is no longer ignored, introducing 
a mesoscopic energy scale AE m = /cBT m . This energy 
scale — in the simple case of longitudinal vibrations — is 
determined by the transverse dimensions a and b of the 
beam through the typical energy difference between suc- 
cessive electronic band minima, given by Eq. (17). It is 
on the order of h 2 ir 2 /2ma 2 to h 2 ir 2 /2mb 2 , to within a 
multiplicative numerical factor arising from the values of 
n y and n z near the Fermi energy. Indeed, Fig. 2 shows 
that when calculated using Eq. (36) which takes 

into account the lateral confinement of the electrons, is 
a function of temperature up to some T m oc 1/a 2 for 
beams with square cross sections. For example, in the 
case of aluminum, the crossover occurs at T m ~ 250k 
for a = b = lOnm (not shown in the Figure), and at 
T m ~ 10k for a = b = 50nm. 

Thus, there are two qualitatively distinct regimes in 
the behavior of electron-phonon damping in our model 
of metallic nanomcchanical resonators. For T sufficiently 
greater than T m , which is obtained either by increasing 
the temperature or by increasing the transverse dimen- 
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FIG. 3. (color online) Schematic plot of the electronic band 
structure near the Fermi energy. The scattering of a longitu- 
dinal phonon with a positive wavenumber q is considered. At 
a very low temperature (T\) the damping is very weak because 
there is no dissipation channel within the narrow temperature 
window around ef- At a somewhat higher temperature (T2) 
the temperature window is wider and the dissipation chan- 
nel with kx, marked by open circles, becomes active. The 
two nearby dissipation channels (marked by open squares) 
remain inactive as they are outside the temperature window. 
The black intervals around the markers represent schemati- 
cally the electronic states that can scatter a phonon when the 
length of the beam is finite and momentum conservation is 
imperfect. 



sions of the beam, the resonator operates in a macro- 
scopic or bulk-like regime, where damping approaches its 
temperature-independent value, as calculated for uncon- 
fined electrons. In the macroscopic regime electrons be- 
have as they do in 3-dimcnsional bulk systems and one 
can safely ignore their confinement by the finite geom- 
etry of the resonator. For T sufficiently less than T m , 
which is obtained cither by lowering the temperature or 
by decreasing the transverse dimensions of the beam, the 
resonator is in a mesoscopic regime, where damping, as 
we explain below, depends non-monotonically on temper- 
ature and other parameters. In this regime the electrons 
behave as a collection of effective 1-dimensional particles, 
each characterized by its own parabolic band, given by 
Eq. (17). 

It follows from Eq. (36), due to the lack of any dis- 
placement in the y and z directions, that a longitudinal 
phonon can be created or annihilated only via intra-band 
electronic transitions. Moreover, due to momentum con- 
servation in the x direction along with energy conserva- 
tion, a phonon of wavenumber q can interact with only 
two specific electronic states in any given band, whose 
momentum in the x-direction is given by k^r, as defined 
in Eq. (37). An electron in the state with k~ can absorb 
a phonon and scatter into the state with fc+, whereas an 
electron with fc+ can emit a phonon and scatter into the 
state with k~ . The difference between these two pro- 
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FIG. 4. (color online) Inverse quality factor of the longitudinal 
mode q = 7r/10 3 nm~ 1 of an infinitely long beam with a = b = 
20nm plotted as a function of temperature, showing that up 
to T ~ 5k electron-phonon damping is dominated by a single 
dissipation channel as explained in the text. 



cesses, for each such pair of electronic states, constitutes 
a dissipation channel, which may contribute to electron- 
phonon damping. It then follows from the form of the 
Fermi-Dirac factors n(l — n), whose width is of the order 
of k^,T, that if the two states happen to lie on either side 
of the Fermi energy, the corresponding dissipation chan- 
nel is active at any temperature. More generically, the 
channel is active only if the two states are located within 
a window of order k^T near the Fermi energy. This is 
demonstrated schematically in Fig. 3. 
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FIG. 5. (color online) Inverse quality factor of the flexural 
mode q = 87r/(0.5 x 10 3 )nm _1 in an infinitely long beam 
with square cross sections a = b = 10, 20, 30, 40, 50nm, and 
rectangular cross sections ai = 40nm, b\ — 20nm and a-2 — 
lOnm, 62 = 20nm, as a function of temperature. Note the 
independence of the asymptotic value on a. 
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In generic situations, at a sufficiently low temperature 
all the dissipation channels are inactive and electron- 
phonon damping vanishes, as can be seen in all the curves 
in Fig. 2. As the temperature increases, at some point 
the dissipation channel that is closest to the Fermi energy 
becomes active, and electron-phonon damping is domi- 
nated by this single dissipation channel. To confirm this, 
we compare in Fig. 4 the total value of Q _1 , as obtained 
from the full sum in Eq. (36), to the contribution of just a 
single term in the sum, corresponding to the dissipation 
channel that is closest to the Fermi energy. The contri- 
bution of a given dissipation channel to the overall damp- 
ing is maximal at a temperature that is roughly given by 



These maxima are 



well pronounced unless a second dissipation channel hap- 
pens to be relatively close to the Fermi energy (as is the 
case for a = 30nm in Fig. 2). As the temperature is fur- 
ther increased the number of active dissipation channels 
gradually increases until eventually the behavior crosses 
over to the bulk regime. 



FIG. 6. (color online) Zoom-in on the low temperature be- 
havior of Q" 1 of the flexural mode given in Fig. 5. 



B. Damping of flexural vibrations in an infinite 
beam as a function of temperature and wavenumber 

The inverse quality factors for flexural modes of a given 
wavenumber q in infinitely long beams — as calculated us- 
ing Eq. (41), taking into account the lateral confinement 
of the electrons — are plotted in Fig. 5 as a function of 
temperature. We again identify two qualitatively dis- 
tinct regimes of behavior — a low-temperature mesoscopic 
regime, where electron-phonon damping is dominated by 
just a single or a few dissipation channels, and a high- 
temperature macroscopic regime, where electron-phonon 
damping approaches its temperature-independent value 
as expected for unconfined bulk electrons. We show be- 
low that the high-temperature limit of the expression for 
Q^ 1 in Eq. (41) is indeed temperature-independent and 
independent of a, and that for long-wavelength modes, 
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with qb <C 1, it behaves approximately as (qb) 2 . This 
is indeed the behavior observed in Fig. 5. The meso- 
scopic regime, which is shown in greater detail in Fig. 6, 
is again characterized by electron-phonon damping that 
tends to zero at very low temperatures, followed by peaks 
at slightly higher temperatures that may result in damp- 
ing that is several times larger than the high-temperature 
constant asymptotic value. 

For flexural modes, it follows from Eq. (41) that the 
emission and absorption of phonons involve inter-band 
electronic transitions, where the value of n z increases or 
decreases by an odd integer. Momentum and energy con- 
servation again restrict the number of electronic states 
that can interact with a given phonon of wavenumbcr 
q, leaving two possible dissipation channels for each pair 
of bands [n y ,n 2 ] and [n y ,n' z ], with n z + n' z an odd in- 
teger. The contribution of these dissipation channels to 
Q _1 decreases rapidly as the difference between n z and 
n' z increases, due to the prefactor n 2 n' 2 / (n 2 . — n' 2 } in 
Eq. (41). Thus, the major contribution to Q^ 1 comes 
from dominant channels with n' z = n z ± 1, while the 
remaining channels can usually be neglected. 

As before, a dissipation channel is active when its two 
electronic states are within a temperature window of the 
order of k^T around the Fermi energy, yet for flexural 
modes the situation is slightly more involved. To see this, 
let us examine the dissipation channels in more detail. 
It follows from Eq. (41) that an electron in the band 
labeled [n y , n z ] with momentum in the x direction given 
by k\ = k + [as defined in Eq. (40)] and a corresponding 
energy e\ [according to Eq. (16)], can absorb a phonon 
with a wavenumber q and scatter into a state in the band 
[n^n^,] with momentum in the x direction k2 — n+ + q 
and energy £2 = £1 + hXu q . The reverse is possible via 
phonon emission. Alternatively, an electron in the band 
labeled [n y , n z ] with momentum in the x direction given 
by &3 = K— and energy denoted by £3, can emit a phonon 
with a wavenumber q and scatter into a state in the band 
[n y ,n' z ] with momentum in the x direction foj = k_ — q 
and energy £4 = £3 — hio q . The reverse occurs via phonon 
absorption. 

Such a pair of dissipation channels, associated with a 
given pair of bands, are depicted schematically in Fig. 7 
for dominant channels with n' z = n z + 1. The energy scale 



Ae(n y 



-associated with the difference in energy be- 



tween the two dissipation channels — is given by 



Ae(n y , \ 



HGtt 2 



(2n z + 1) oc 



1 



(43) 



whereas the average energy e(n y ,n z ) associated with 
both channels, in the limit of long- wavelength vibrations, 
is given by 



e(n y , n z ) = lim e. 



<7&->0 



H 2 7T 2 

2m 



JT 

'qb 2 



(44) 

where Ei is any one of the four energies £\ ... £4 mentioned 
earlier and shown schematically in Fig. 7. 
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FIG. 7. (color online) A schematic plot of two electronic en- 
ergy bands characterized by the same n y and with n' z = n z + 1 
showing the pair of dissipation channels associated with a flex- 
ural mode with wavenumber q > as explained in the text. 
At the temperature T, indicated by a dashed line, only one 
of the two channels is active; at a higher temperature both 
channels would be active and one would lose the ability to 
resolve them, as demonstrated below in Fig. 10. The black 
intervals around the four markers represent schematically the 
electronic states that can scatter a phonon in finite-length 
beams where there is only approximate conservation of mo- 
mentum. 



At high temperatures, the leading term in an expansion 
of Eq. (41) in powers of f3ftw q <C 1, shows that each active 
pair of dominant dissipation channels with n' z = n z ± 1 
is responsible for a contribution to the damping that is 
proportional to f3huj q /abq, which in turn is proportional 
to flq/a. To estimate the total damping, one simply 
needs to count the number of integer pairs (n y ,n z ) cor- 
responding to dissipation channels whose electronic ener- 
gies s(n y , n z ) are a distance of order fcsT within £f- The 
density of such states, obtained directly from Eq. (44), 
yields a number that is proportional to qab 2 k^T, and 
therefore to a total damping proportional to {qb) 2 and 
independent of a and T, as mentioned earlier. 

In Fig. 8 we plot Q _1 as a function of phonon wave- 
number q for several temperatures. At low temperatures, 
in the mesoscopic regime, Q~ l fluctuates greatly as q is 
increased, although the general trend is an increase in 
the damping. However, as the temperature crosses over 
to the macroscopic regime the fluctuations are suppressed 
and the damping converges to a single monotonically ris- 
ing curve. The large fluctuations in Q^ 1 in the meso- 
scopic regime again stem from the small number of active 
dissipation channels, whose positions depend very sensi- 
tively on the value of q. Thus, for narrow temperature 
windows around the Fermi energy, small changes in q 
can easily cause an active dissipation channel to become 
inactive or vice versa. On the other hand, in the macro- 
scopic regime the number of active channels is large, and 
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FIG. 8. (color online) Inverse quality factor of flexural modes 
with q = 71-/(0.5 x 10 3 )nm _1 to 60 x 71-/(0.5 x 10 3 )nm _1 (which 
approximately correspond to the first 60 flexural modes of a 
I = 0.5/^m beam) in an infinitely long beam with a — b — 
20nm for different temperatures. 

small changes in q affect the damping only monotonically. 
Fig. 8 confirms that for small q in the macroscopic regime 
(J -1 is indeed proportional to q 2 . 



C. Damping of flexural vibrations in an infinite 
beam as a function of ef 

As an example of the delicate sensitivity of electron- 
phonon damping to changing system parameters in the 
mesoscopic regime, we consider the dependence on the 
Fermi energy. The latter may be varied in real experi- 
ments by means of an external gate voltage, suggesting 
an intriguing way to probe the electronic spectrum by 
measuring the Q of mechanical vibrations. As for lon- 
gitudinal modes, the existence of individual dissipation 
channels can be discerned as long as k^T is smaller than 
the energy scales characterizing the difference between 
channels. For flexural modes there are two such energy 
scales. The first is the energy difference (43) between two 
dissipation channels belonging to the same pair of energy 
bands; and the second and larger energy scale is the dif- 
ference between pairs of dissipation channels belonging 
to different pairs of electronic energy bands (44) . 

Low-temperature results for T = 10k are shown in 
Fig. 9 for flexural vibrations in an infinite beam with a 
square cross-section. The damping is zero up to a thresh- 
old value of £p, which according to Eq. (44) is approxi- 
mately equal to e(l, 1) ~ 97r 4 ?i 2 /8mq 2 & 4 . As the Fermi 
energy is raised above this threshold, the damping ex- 
hibits a series of narrow peaks, while remaining zero for 
most values of £v ■ For the particular beam geometry and 
energy range considered in Fig. 9, the observed damping 
peaks arc associated with pairs of dissipation channels 
belonging to bands [n y , n z = 1] and [n y , n' z = 2], with n y 
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FIG. 9. (color online) Inverse quality factor of the flexural 
mode q = 27r/10 3 nm _1 in an infinitely long beam with a = 
b = lOnm at T — 10k as a function of the Fermi energy. 
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FIG. 10. (color online) Inverse quality factor of the flexural 
mode q = 27r/(10 3 )nm _1 in an infinitely long beam with a = 
b = lOnm at T — 2k and 10k as a function of the Fermi 
energy. 



increasing by 1 from one peak to the next. The increase 
of n y by 1 between consecutive peaks results in a linear 
increase of the energy separation between the peaks as 
the Fermi energy increases, in accordance with Eq. (44). 

In Fig. 10 we zoom-in on one of the peaks and lower 
the temperature even further to reveal the expected in- 
ner structure, consisting of two sub-peaks. The two sub- 
peaks are distinct from each other if the temperature is 
small compared to the energy difference Ae between the 
two energy channels as is the case for T = 2k in Fig. 10. 
As the temperature increases the sub-peaks become less 
distinct — as can be seen in Fig. 10 for T = 10k — and 
eventually merge as the temperature becomes compara- 
ble to the energy difference between the adjacent chan- 
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FIG. 11. (color online) Inverse quality factor of the flexural 
mode q = 27r/10 3 nm _1 in an infinitely long beam with a = 
b = lOnm at T = 300k as a function of the Fermi energy. 

nels as given by Eq. (43). 

In Fig. 11 we show how Q _1 varies with £p at a 
high temperature of T = 300k. The dense peaks near 
tF = 14eV in Fig. 9 broaden and merge into a single large 
asymmetric peak. As the Fermi energy increases — and 
along with it so does the separation between the probed 
pairs of dissipation channels — individual peaks start to 
emerge and become visible. However, the sub-peak struc- 
ture within each main peak remains indiscernible. 

Eqs. (43) and (44) continue to provide a good qualita- 
tive understanding of the calculated damping peaks for 
other system parameters, as long as qb remains smaller 
than 1. For example, doubling the transverse dimensions 
of the beam to a = b = 20nm lowers the threshold en- 
ergy by a factor of 2 4 = 16 to somewhat below leV, and 
the entire peak structure shown in Fig. 9 is shifted to 
between about leV and 2.5eV. At higher Fermi energies 
the structure becomes less ordered than the one shown 
in Fig. 9, because the energy is sufficiently large for the 
value of n z to start changing between peaks as well. 

D. Electron-phonon damping in the presence of 
other dissipation mechanisms 

The calculation of the phonon decay rate in Eq. (24) 
assumes that initially both electron and phonon life- 
times are infinite, and therefore their spectral functions 
are given by Dirac delta functions. This assumes that 
the electron mean-free path is greater than the wave- 
length of the vibration mode, and that all other damping 
mechanisms are much weaker than the damping caused 
by electron-phonon scattering. As a quick aside, we 
wish to account for the effect of other damping mech- 
anisms by adapting our calculation to allow for phonon 
states with a priori finite lifetimes, or broadened spectral 
functions. We adopt an approximate phenomenological 
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FIG. 12. (color online) Variation of Q _1 of the flexural mode 
q = 27r/10 3 nm~ 1 in an infinite beam with a = b = lOnm at 
T — 10k as a function of the Fermi energy. A solid (blue) line 
shows the results obtained if one ignores other damping mech- 
anisms, while a dashed (red) line shows the results obtained 
using a broadened phonon spectral function in the form of a 
Lorentzian with Qq" 1 = F/hug = 10 -3 . 

method used in electron transport simulations, 8 '~ 90 in 
which collisional broadening of electrons is taken into ac- 
count. We replace the sharp energy conservation delta 
functions in the expression for the decay rate (24) with 
a convolution over a broadened spectral function of the 
phonon in the form of a Lorentzian, 

A(s)= 1 - hT ' 2 2 . (45) 

We ignore any shift in the energy of the phonon that 
might be caused by other damping mechanisms and 
continue to take the electron spectral functions to be 
delta functions. Thus, we replace the delta functions 
5 (£k — £k' i Hujq) in Eq. (24) with their convolution with 
the Lorentzian spectral function (45) of the phonons, 

A (e k - £k' ± tiojq) 

= J A (E - huj q ) 5 (e k - e k < ± E) dE. (46) 

We assume an initial broadening with Q^ 1 = T jhijj q 
that is much larger than the values calculated above for 
electron-phonon damping. This corresponds to situa- 
tions where other sources of damping completely obscure 
the effect of electron-phonon damping. Nevertheless, if 
for such initial broadening the contribution of electron- 
phonon damping is not significantly altered, then it will 
not be affected when the initial damping is smaller and 
electron-phonon damping is a dominant effect. As an ex- 
ample, we take Qq 1 to be equal to 10~ 3 (which reflects 
a metallic beam with a relatively poor quality factor) 
and recalculate one of the damping peaks obtained in 
the previous section as function of Ef for a flexural mode 
with q = 27r/10 3 nm~ 1 in an infinitely long beam with 
a = b = lOnm at T = 10k. As can be seen from Fig. 9, 
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FIG. 13. (color online) Inverse quality factor of the lowest 
forty flexural modes of a finite beam with I — 0.5/^m and 
a = b = 20nm at T = 10k. The results obtained for the same 
values of q for an infinite beam are given for comparison. 



Q q 1 of this mode is smaller by at least an order of magni- 
tude compared to Qq 1 - In Fig. 12 we show that the effect 
of an initial broadening on such a peak is very small — a 
result which holds also for the rest of the peaks as well 
as for low temperatures. This indicates that even for 
relatively high damping by other dissipation mechanisms 
(compared to "stand-alone" electron-phonon damping), 
the contribution of electron-phonon interaction to the to- 
tal damping is nearly unaffected. 



E. Damping of flexural vibrations in a finite beam 

Finally, we consider the effect of the finite length of the 
beam by calculating the inverse quality factors of flexural 
modes in a finite beam using Eq. (42), which takes into 
account the lateral confinement of the electrons. Fig. 13 
shows the low-temperature behavior of Q^ 1 as a function 
of wavenumbcr q in a finite beam, as compared with the 
values calculated earlier for an infinite beam. We can 
still identify the non-monotonic dependence on parame- 
ters, which is characteristic of the low-temperature meso- 
scopic regime, although the sharp features average out in 
the finite beam. The damping in the finite beam devi- 
ates significantly from the damping in the infinite beam, 
exhibiting smaller fluctuations away from the general in- 
creasing trend, although the relative difference between 
the two results decreases with increasing wavenumber q. 

The difference between the finite beam and the infinite 
beam is due to the replacement of the momentum delta 
functions in Eq. (20) for an infinite beam with the C n 
functions appearing in Eq. (22) for a finite beam. Thus, 
exact momentum conservation is replaced by an approx- 
imate or an imperfect momentum conservation. As a 
consequence, the discrete electronic states that form a 
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FIG. 14. (color online) Inverse quality factors of the lowest 
forty flexural modes of the same beam as in Fig. 13, but at 
T = 300k, compared again with the inverse quality factors of 
an infinite beam. 

dissipation channel in an infinite beam are replaced by 
small intervals of states — shown schematically in Fig. 7 — 
causing the smearing or averaging effect that we observe. 

In order to recover the infinite beam results, the C n 
functions need to be sufficiently sharp compared to all 
other scales related to the momentum conservation in 
the x direction. At high temperatures this is simply the 
case for short-wavelength modes, where q n l ^> 1. How- 
ever, at low temperatures we need to require in addition 
that the C n functions be sufficiently narrow compared to 
the width of the Fermi-Dirac factors of the form n(l — n). 
Thus, as the temperature increases and with it the width 
of the Fermi-Dirac factors, the difference between the 
damping of the finite beam and of the infinite beam de- 
creases, as shown in Fig. 14. 

VI. SUMMARY AND CONCLUSIONS 

We have estimated the effect of free electrons in ther- 
mal equilibrium on the damping of mechanical vibra- 
tions in metallic nanoresonators in the form of thin long 
beams. In doing so we made a sequence of assumptions: 
(1) phonon wave functions were obtained by quantizing 
the normal modes of vibration obtained by the simplest 
thin-beam elasticity theories; (2) electrons were consid- 
ered as free noninteracting particles either totally un- 
confined, as in bulk material, or confined by the lateral 
surfaces of the resonators; (3) electron-phonon interac- 
tion was described by means of a short-range electro- 
static potential, where electrons couple to phonons only 
through local volume changes, which are largest for lon- 
gitudinal modes leading to Q's on the order of 10, are 
smaller for flexural modes leading to Q's on the order of 
10 3 to 10 6 , and are non-existent for twist modes; (4) the 
calculation itself was performed at the level of Fermi's 
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Golden Rule, considering an electron-phonon system at 
thermal equilibrium with just a single additional phonon 
in the vibration mode of interest; (5) the a priori lifetimes 
of the electrons and the phonons were assumed infinite, 
thereby imposing the adiabatic or high-frequency limit, 
Lo q T e > 1, and assuming that all other vibration damping 
mechanisms are much weaker (although we showed that 
the latter requirement is not too stringent); and (6) to ob- 
tain actual results we used the material properties of bulk 
aluminum, neglecting their dependence on size and tem- 
perature. Aluminum was chosen because it is expected to 
exhibit a relatively high degree of electron-phonon damp- 
ing among the common metals used to fabricate nanomc- 
chanical resonators. Despite all of our assumptions and 
approximations, we believe that our results provide an 
adequate qualitative understanding of the behavior of 
electron-phonon damping in metallic nanorcsonators in 
the adiabatic limit. 

Within these assumptions, electron-phonon interaction 
is restricted by energy conservation (through Fermi's 
golden rule), and by momentum conservation along the 
length of the beam when the beam is considered to be 
infinite. These conservation laws severely restrict the 
number of electronic states that can contribute to the 
damping of a given vibration mode, leading to the ap- 
pearance of discrete dissipation channels. The energy 
separation between dissipation channels is determined by 
the dimensions of the cross-section of the beam, and for 
flexural modes also by the phonon wavenumber. This in 
turn sets a temperature scale T m for the crossover from 
high-temperature macroscopic behavior, where electron- 
phonon damping behaves as if the electrons were in the 
bulk, to low-temperature mesoscopic behavior, where 
damping is dominated by just a few dissipation channels 
and exhibits sharp non-monotonic changes as parame- 
ters are varied. As a consequence, for example, at suf- 
ficiently low temperatures one can perform spectroscopy 
of the electronic states by measuring the quality factor 
of a given vibration mode, while varying an external gate 
potential. 

In beams of finite length one lacks translation invari- 



ance, and momentum conservation along the length of 
the beam is lost. Nevertheless, for sufficiently long beams 
compared to the wavelength of the vibration, one still 
retains approximate, or imperfect, momentum conser- 
vation, where delta functions are replaced by finite yet 
sharply peaked functions, whose width is inversely pro- 
portional to the length of the beam. The infinite-length 
limit is recovered when the width of these functions be- 
comes sufficiently small compared to all other relevant 
momentum scales along the length of the beam. At high 
temperatures the only relevant scale is the wavelength of 
the mode, and the infinite-length limit is simply recov- 
ered when the wavelength becomes small compared to 
the length of the beam. At low temperatures the infinite 
limit is obtained only when, in addition, the width of the 
imperfect conservation functions becomes smaller than 
the widths of the Fcrmi-Dirac factors that appear in the 
expressions for the inverse quality factor. Thus, longer 
beams are required as temperature is decreased. 

To conclude, we have demonstrated several interesting 
features of the damping of vibrational modes in metal- 
lic nanorcsonators. These features arise from finite size 
effects that lead to a quantization of the electronic spec- 
trum as well as to imperfect momentum conservation, 
and are dominant in the low-temperature mesoscopic 
regime. At high temperatures the damping becomes sim- 
ilar to the one obtained for unconfined electrons. Our 
analysis was carried out within a simple specific geome- 
try. Nevertheless, we believe that the different regimes we 
have identified are relevant to electron-phonon damping 
in other types and geometries of metallic nanorcsonators. 
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